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Abstract 

We conjecture exact and simple formulas for some physical quantities in two quantum chains. 
A classic result of this type is Onsager, Kaufman and Yang's formula for the spontaneous mag- 
netization in the Ising model, subsequently generalized to the chiral Potts models. We conjecture 
that analogous results occur in the XYZ chain when the couplings obey Jy + Jy Jr + JxJz = 0, 
and in a related fermion chain with strong interactions and supersymmetry. We find exact for- 
mulas for the magnetization and gap in the former, and the staggered density in the latter, by 
exploiting the fact that certain quantities are independent of finite-size effects. 

Onsager's computation of the exact partition function of the two-dimensional Ising model HI 
is one of the great triumphs of theoretical physics. This result now can be reproduced easily, by 
using Kaufman's mapping of the spins to free fermions Q. The computation of the spontaneous 
magnetization, by Onsager and Kaufman f3l and by Yang ||4l, is a second triumph: because the map 
from spins to fermions is non-local, the computation was and remains quite intricate |j5||. Their final 
result is exceptionally simple. The spontaneous magnetization in the ordered phase k < lis exactly 
(1 - A;2)V8 in the large-lattice limit; 1/k = sinh(2 J/fe^T) sinh(2 J7A;Br), where J and J' are the 
usual Ising couplings for the horizontal and vertical links of the square lattice. 

It is natural to guess that the simplicity of this formula is a consequence of the model's underlying 
free-fermion nature. Thus it is remarkable that an elegant generalization of Onsager, Kaufman and 
Yang's formula occurs in a series of models most decidedly not free fermions. The chiral Potts model 
is a parity -breaking Z^r generalization of the Ising model with some amazing properties ||6l|7l. One 
is that the order parameters for spontaneously breaking the 7,^ symmetry are given by a formula just 
like the Ising model, as conjectured in ref. [8] and proved more than 15 years later in a tour deforce 
of Baxter's ||9l. Labeling the spin at site j by a variable aj = . . . N — 1, the exact result as the 
number of sites goes to infinity is 

^^2-,Tiraj/N'^ = (1 _ ;,2^)r(Ar-r)/(2Ar2) 

The lattice parameter k in ([T]) is not renormalized: it is a coefficient of one of the terms in the 
corresponding quantum Hamiltonian. Nevertheless, the expression for the order parameters in (O 
is exact for any value of k, ranging from the critical point A; = 1 to the completely ordered point 
k = 0. This is unusual even for integrable models; when order parameters can be computed they are 
typically given by elaborate combinations of elliptic theta functions (see ref. LIOJ ). 

In this paper we conjecture exact formulas analogous to ([U in two quantum chains with strong in- 
teractions: the XYZ chain along a special line of couplings |[m[T2l[T3l[T4l . and interacting fermions 
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with supersymmetry lfT5l[T6]| . The conjectures result from studying series expansions around a triv- 
ially solvable limit, the analog of A: = above. We find that for a system with L sites, the terms in 
these expansions up to order L are independent of L. We refer to such quantities as scale free. We 
can thus compute them exactly by finding the ground states explicitly for small systems. The analogs 
of ([Hi then are obtained by summing the series. Since this yields the correct critical exponents for the 
model, this provides strong evidence that the conjecture is exact in the L — > oo limit. 

To motivate our study, we note two special properties of the chiral Potts model. One is that along 
a line in parameter space, it possesses a useful symmetry algebra, known as the Onsager algebra 
im 111, which allows the explicit construction of an infinite sequence of conserved quantities. A 
second (under-appreciated) property is that in the corresponding field theory in the scaling limit, the 
coefficient of the Lorentz-symmetry breaking perturbation does not renormalize ifTTl . 

Supersymmetric field theories also possess such special properties. Because the Hamiltonian 
is part of the supersymmetry algebra, supersymmetry does much more than just grouping of states 
into multiplets. One can often prove the existence of zero-energy ground states by computing the 
Witten index lITSl . Moreover, in some cases there are non-renonnalization theorems. For example, in 
the scaling limit of the models described below, the superpotential does not receive any corrections 
beyond tree level in perturbation theory |[T9l . This means that some physical quantities (for example, 
the gaps of certain kink states) depend simply on the parameters in the Hamiltonian. 

This motivates us to study quantum chains whose scaling limits aix described by supersymmetric 
field theories. Our first example is a special case of the well-known XYZ chain ifTOl . The Hilbert 
space (C^)®^ is a two-state system at each site on the chain, and the Hamiltonian is 



where the o"" are the Pauli matrices and Eq is a constant. For now we take periodic boundary con- 
ditions, so that o-£^^ = erf. When Jx = Jy, the Hamiltonian preserves the numbers of up spins and 
down spins individually; elsewhere these numbers ai^e only conserved mod 2. For L odd, all states 
including the ground state are therefore paired by flipping all the spins. If one of the Ja vanishes, 
the chain can be mapped onto free fermions by the usual Jordan- Wigner transformation; otherwise, 
the mapping gives interacting fermions. Whenever Jx = Jy and \Jz\ < Jx (and values related by 
permuting the Ja), the model is critical, and is called the XXZ chain. Along this critical line, a 
free-boson field theory describes the scaling limit. Neai^ this critical line, it can be described by the 
sine-Gordon field theory. 

The field theory of the XYZ chain is supersymmetric along a particular line in its two-parameter 
space (see e.g. ref. Il20l ). Because the chain is integrable ITOl . it is easy to identify the supersymmetric 
critical point in the XXZ chain: it is at = —\Jx\/2. The XXZ chain here has many fascinating 
properties (see e.g. ref. Il2l1 ). In fact, long ago Baxter found a simple formula for the exact ground- 
state energy as L — oo along the entire line 



Namely, for £^0 = —{Jx + Jy + Jz), the ground-state energy along this Une goes to zero as L — 00. 
Moreover, it was conjectured that the lowest eigenvalue of Hxyz along this line is exactly zero when 
L is odd 11211 . just as in supersymmetric models. This was subsequently proved (for L odd as well 
as for L even with twisted boundary conditions) at the critical point by showing the XXZ chain has a 
hidden supersymmetry relating chains with different numbers of sites |[T6ll22l . Moreover, there are a 
host of other fascinating and special results along this line l[T2l[T3l[T4l . all reminiscent of the special 
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results occurring in fermion chains with an explicit supersymmetry lfT5l[T6ll23]| . Thus the XYZ chain 
along the line Q indeed should correspond to a supersymmetric field theory in the scaling limit; for 
this reason we dub this the sXYZ chain. 

Our second chain is a staggered version of a fermion chain with a built-in supersymmetry lfT5l[T6ll . 
These models are defined from the supersymmetry operator Q obeying = 0. The Hamiltonian 
H = QQ'^ + Q"^ Q the commutes with Q and We study the supersymmetric Hamiltonian acting 
on a Hilbert space spanned by spinless fermions, with the additional restriction that fermions may not 
be on adjacent sites. The supersymmetry operator in terms of fermion creation operators cj is 

Q = ^ Aj(l - rij-i)(l - ?ij+i)cj, (4) 
i 

where n-j = c^jCj. Q squares to for any choice of the complex numbers Xj, so the Hamiltonian 



3/ 

-ffssF = X] I (1 - nj_i)(A*Aj+ic]cj+i + /i.c.)(l - nj+2) + |Aj|^(l - nj_i)(l - n^+i) 
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is supersymmetric. The first term allows hopping preserving the restriction, and the second is com- 
prised of a chemical potential and a next-nearest-neighbor repulsion. When the number of fermions 
is / and the number of sites is 3/, for periodic boundary conditions there are two ground states for 
any values of the Xj ifTSl . Here we consider the staggering X^i = Asj+i = 1 and A3j+2 = z, and 
so we label this model ssF (for supersymmetric staggered Fermions). The Bethe equations for the 
unstaggered case z = 1 and for the critical sXXZ chain are the same up to boundary conditions |[T6l . 
and so the critical field theories must be the same. Staggering the model perturbs it away from this 
critical point, and since there is only one Lorentz-invariant supersymmetry-preserving perturbation, 
its scaling limit should be the same as the supersymmetric field theory describing sXYZ. 

For the remainder of this paper we describe some of the remarkable properties of these models. 
The key to much of our analysis is to expand various quantities around a limit where the model can 
be solved trivially. An amazing property of these chains is that for certain quantities, the coefficients 
of the terms in this expansion are scale free. 

We parametrize the sXYZ line dS]) by 

Ja: = 2s{s-3), Jy = 2s{s + 3), = 9 - s^, (6) 

so that £"0 = 3(s^ + 3). The critical points are at s = ±1, oo, while at the trivially solvable points 
s = 0, ±3 only one of the three terms in Q remains. At s = 0, only / 0, so the ground 
state |0) has all spins the same. The spontaneous magnetization per site Ml{s) = (0|(j||0), obeys 
Ml{0) = 1 in the sector with an even number of down spins. We find, by using Maple to compute 
the exact ground state, that the power-series expansions of the magnetization for odd L are 

Mg = l-4s2-12s^ + 188s6-844s*^ + 380si° + ... 

M7 = l-4s2-12s^-52s^ + 2516s*^- 18004s^° + ... 

Mg = l-4s^-12s^-52s^- 284s^ + 33516s^° + ... 

Mil = l-4s2-12s^-52s^- 284s^- 1764si° + ... 

where s = s/3. The trend is obvious: the order terms in the expansion are independent of L when 
n < L. The magnetization appears to be scale free near s = 0. Doing this to L = 17 yields what is 
presumably the exact expansion as L — 00: 

Ml{s) = 1- 4s^ -12s^ -52s^ - 284 - 1764 - 118205^^ - 832205^'' - 606780s^'' + . . . 
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To understand how to sum the series and find a simple formula Moo{s), we examine the expected 
behavior at the critical point s=l. The dimension of the "thermal" operator that perturbs away from 
s=l onto the sXYZ line is 4/3, while the dimension of the magnetization operator is expected to be 
1/3 Uni. Indeed, the finite-size values at criticality fitnicely to Mz,(l) w .95527 L-^/^{l +0(1'^ )). 
Thus as s 1", Moo(s) should vanish as (1 - s)^ with /3 = (l/3)/(2 - 4/3) = 1/2. This square- 
root singularity suggests looking at the series expansion of {Ml{s))^: 

{ML{s)f = 1 - - - Bs^ - - 8^1° - 85^2 _ g^-M _ g-ie ^ q^^l+I) 
Summing this series yields a conjecture for the exact magnetization in the ordered phase s < 1: 

^-(«) = 3 (^3^J • (7) 

We emphasize that we do not assume anything about the behavior at the critical point s = 1; the 
only role of the scaling argument is to suggest that we square M. The fact that the expected critical 
behavior for the magnetization emerges from the expansion around s = is to us a compelling 
argument that the formula ([7]) is exact. 
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Figure 1: Ml{s) for s < 1; the solid red curve is the conjecture for Moo(s), while the dashed curves 
are for L = 5, 9, 13, 17. 



We plot this function and the finite-size curves in figure [T] Even with the finite-size effects 
near s = 1, it is clear that the finite-L curves are approaching the conjectured curve. A numerical 
calculation using the iTEBD method gives the same curve to high accuracy f24|. Moreover, with a 
change of variables, O gives the same "gr-series" obtained by exploiting the integrability of the chain. 
The corresponding quantity in the eight- vertex model, the spontaneous polarization, is |[25ll26l 



4vr / ^.(0,gV(^-^'^/-)) y 
° 2^-^\ ^2(0,<?V2) j W 

when written in terms of Jacobi theta functions |[27l . The sXYZ line ^ corresponds to setting the 
crossing parameter rj = 7r/3. Using the Boltzmann weights of the eight- vertex model at this point. 
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we find ^ 

for < s < 1. Inserting this into ^ yields the same expansion in q as that of ([S]) to ~ 500 terms. 
We have found other exact formulas using these methods. Letting Hj = a'^a'^j^^, for s < 1 

(0|//||0) = l + Af{-l + s^ + ?>s'^ + bs^ + 7s^ + ...) 

We can find this expectation value for s > 1 by expanding around s = 3, where Jx = Jz = 0. 
Letting t = (3 - s) /6, 

{0\H^\0) = ^{2t + 3t'^ + 4t^ + 5t^ + ...) 

Using the g-series representation ^ we showed analytically that (fTOl) matches the exact results |[28i . 
The expectation values of Hj and Hj can be found from these by using the duality symmetries 
s — )• (3 — + 1) and s — )• —s, or by using Hellmann-Feynman theorem (0|dHsXYz/ds|0) = 0. 

Not only ground-state properties are scale free: the gap is as well, and obeys an elementary 
formula. To define the gap, we exploit the fact that there is a spontaneously broken Z2 symmetry 
away from the critical points. It is thus natural to think of the gapped excited states as kinks sepai^ating 
regions of the two ground states. This picture is supported by the computation of the exact scattering 
matrix for these kinks in the supersymmetric field theory |[20l . For an odd number of sites and 
periodic boundary conditions, we expect the lowest-energy excited states to be two-kink states. Since 
the kinks interact, the energy is less than twice the kink gap. Thus to define the gap to the one- 
kink state, we consider an even number of sites with twisted boundary conditions (a spin-flip defect): 
(T^_l_^ = —erf, cr^_^-^ = —CTi, and cr,^^i = +(Tf' . Near s = 0, the interactions away from the boundary 
favor aUgning the spins, but the twist forces the energy to be order J^. 

The result for the gap is more transparent when we rescale the Hamiltonian H — )• H/s^, and 
consider the region between the trivially solvable point s = 3 and the critical point at s — )• 00; the 
gap in other regions is obtained by exploiting duality. We found the exact one-kink energy A for 
sizes up to L = 10. Expanding this in a power series around s = 3 in terms of f = 1 — 3/s, we find 

Al = 4-67; + 3t;V2 + fV4 + 3w^/32 + ... 

= 4 0y^' + O(t;^/2) (12) 

Thus at the critical point s — 00, the gap vanishes with exponent v = 3/2. This is exactly what one 
expects with dimension-4/3 thermal operator: v = 1/(2 — 4/3) = 3/2. 

We now turn to the supersymmetric staggered fermion model with Hamiltonian and show that 
not only does it possess scale-free quantities similar to those of the sXYZ chain, but that the models 
are deeply related on the lattice, not just in their scaling limits. For / fermions on 3/ sites, H^sf has 
two zero-energy ground states like HsXYZ- Here, however, the two ground states ai^e not related by 
symmetry: because the fermions cannot occupy adjacent sites, there is no analog of spin-flip symme- 
try. To define basis vectors for the two-dimensional space of zero-energy states unambiguously, we 
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exploit the parity symmetry j — )• 3/ + 1 — j. One ground state, denoted is even under parity, 
while the other ground state |— )z is odd. These two ground states are quite different from each other, 
as is easy to see by studying them in the solvable limits z — 0, oo. Since H = QQ^ + Q^Q, any 
zero-energy ground state must be annihilated by both Q and Q^. Letting | j) label the three states with 
a fermion on every third site 3i + j, we have at z — oo, |±)oo = ± |3))/\/2. In the other limit, 
|+)o = 1 2), but the odd-parity ground state is a sum over all configurations without a fermion on the 
sites 3i + 2: |-)o = nf=i(4+i " 4)|empty)/2//2 . 
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Figure 2: D~^{z) — D (z); the solid red curve is conjecture (fT4l) . and the dashed curves are for sizes 
3/ = 12,18,24. 

We find exact formulas for the staggered fermion densities D^{z) = z{i^\c\i_iC3i-^i\±) z- These 
have been studied at the critical point z = I, using numerics |[23l and using conformal field theory 
1291 . and we extend these results to all z. In the solvable limits, we have D^'(oo) = D~{oo) = 
D~{0) = and D^{0) = 1. Moreover, at the critical point z = 1, the full translation symmetry of 
the model is restored. This requires that D^{1) + D~{1) = 2/3. We thus expect that — D~ 
behaves like the magnetization in the sXYZ chain, vanishing as / — cxo when z > 1, but non-zero 
for z < 1. By finding the exact ground state in sizes up to / = 8 (24 sites), we obtain for z small 

+ D~ = 1 - 35^ + 3z'^ -3z^ + 3z^ - ... 

^~'^'\o{z'f), (13) 



8 + z2 

-D- = l-hz^ - 3i* - 29 f - 1315 
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^ 2 +Oiz'-f) (14) 

where z = z/\/^. We see the same square-root singularity in Z)+ — that we did for the magneti- 
zation in sXYZ. The series expansion around z = oo gives 

_^ 2 6 26 134 762 4614 

D +D - ^2 -^ + ^-^ + ^-^ + --- 

^ +0{z~^^) (15) 



z2 + zy/8 + ^2 + 2 
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The finite-size curves for + D~ are almost indistinguisliable from tlie asymptotic form, because 
of tlie exact result at z = 1 and the scale-free behavior in the small- and large-z limits. The curves 
for Z)+ — D~ are plotted in figure |2l 

Field-theory dualities need not be exact in the corresponding lattice models, or can be very subtle 
(e.g. the Kramers-Wannier duality of the Ising model). Thus even though the sXYZ chain has the 
duality symmetry s — )• (3 — s)/(s + 1), the coiTcsponding duality is not obvious in the ssF chain. 
Nevertheless, we have non-trivial evidence that there is such a duality exchanging the \z\ > 1 and 
\z\ < 1 phases. This becomes apparent when we simplify (fTSl ) by defining the new coupling S = 
3z/\/i2 + 8, so that D+ + D' = 2(3 - 5)/(3S' + 3) asymptotically for z > 1. The z = 1 critical 
point is at 5 = 1, while the solvable points 2 = and z = oo con^espond to 5 = and 5 = 3 
respectively, the same as the value of s in the sXYZ chain. We find that 

Z)+ (5) D~{S) = D+ (S) (S) 

for S = {3 — S)/{S + 1). This relation holds for all finite sizes up to 3/ = 24 sites, and of course 
for the asymptotic formulas as well. It thus seems very likely that this is a general symmetry of the 
ssF chain, but we have not yet found the corresponding symmetry of the Hamiltonian. 

The relation between the sXYZ and ssF chains goes even deeper. In a remarkable series of papers 
|[T3l . Bazhanov and Mangazeev showed that (at least for small systems) the ground states themselves 
are related to the tau functions of the Hamiltonian hierarchy of the Painleve VI non-linear differential 
equation. They find a recursion relation for the coefficient of the state with all spins down in the 
wavefunction, normalized so that it is a polynomial in s. This same polynomial appears in the ground 
state in our ssF chain! It appears (up to a convention-dependent overall power of S) for example as 
the coefficient in of the state |2) defined above, when z is rewritten in terms of S. Moreover, the 
normalizations of the ground states are related to the same polynomials, just as in ref. (W]- Thus the 
ssF ground states can be related in the same fashion to Painleve VI. 

We have presented conjectures for exact results in two interacting chains. These include simple 
formulas for the spontaneous magnetization and the gap in the XYZ chain when the scaling limit 
is a supersymmetric field theory. We believe that the evidence for these conjectures is convincing. 
Moreover, these chains quite obviously have a great deal of symmetry structure left to be uncovered. 
In particular, all the evidence - the scale free property, the important role of supersymmetry, and the 
precise relations between the ground states of the two chains - makes it seem likely to us that there 
is an infinite-dimensional symmetry in both models similar to the Onsager algebra lU. Each model 
will then correspond to a different presentation of this algebra. 

This work was supported by the NSF grant DMR/MSPA-0704666. 
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